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Abstract. We demonstrate that the inclusion of the helicity flux in the magnetic helicity balance in the nonlinear 
stage of galactic dynamo action results in a radical change in the magnetic field dynamics. The equilibrium 
value of the large-scale magnetic field is then approximately the equipartition level. This is in contrast to the 
situation without the flux of helicity, when the magnetic helicity is conserved locally, which leads to substantially 
subequipartition values for the equilibrium large-scale magnetic field. 
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1. Introduction 

The large-scale magnetic fields of galaxies are thought to 
be generated by a galactic dynamo due to the simulta- 
neous action of the helicity of interstellar turbulence and 
differential rotation (see, e.g., Ruzmaikin et al. 1988). The 
kinematic stage of the galactic dynamo, i.e. the evolution 
of a weak magnetic field with negligible influence on the 
turbulent flows, seems to be clear, while the nonlinear 
stage of dynamo evolution is a topic of intensive discus- 
sions (for reviews, see Beck et al. 1996, Kulsrud 1999). The 
most contentious issue is the question of the equilibrium 
magnetic field strength at which dynamo action saturates. 

A naive viewpoint is that the saturation level for the 
large-scale magnetic field is given by the equipartition be- 
tween kinetic energy and the energy of the large-scale 
magnetic field B (see, e.g., Zeldovich et al. 1983). The 
motivation is that the equations describing large-scale dy- 
namo action contain the mean, but not the total, magnetic 
field. This naive outlook leads to models of dynamo gen- 
erated magnetic fields which are in basic agreement with 
the available observational information. 

Vainshtein and Cattaneo (1992) formulated a more 
sophisticated argument, suggesting that the equilibrium 
magnetic field should be determined by a balance between 
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the kinetic energy and the energy of the total magnetic 
field. The simplest models of dynamo generation then re- 
sult in the estimate b/B ~ Rm 1 / 2 , where b is the small- 
scale magnetic field, and the magnetic Reynolds number 
Rm rs 10 s for the interstellar turbulence (or even much 
larger if a microscopic diffusivity instead of ambipolar dif- 
fusion is used; cf. Brandenburg & Zweibel, 1995). Thus 
the ideas of Vainshtein and Cattaneo lead to the conclu- 
sion that a dynamo generated large-scale galactic mag- 
netic field must be negligible in comparison with that ob- 
served, and so the generation of the observed field must 
be connected with another mechanism. However, no other 
general and realistic mechanism for galactic magnetic field 
generation is currently available. 

The arguments of Vainshtein and Cattaneo do not 
seem inevitable. For example, a dynamo generated mag- 
netic field can itself produce helicity, so the nonlinear ef- 
fects can even amplify rather than suppress field gener- 
ation at the initial stages of nonlinear evolution (Parker 
1992, Moss et al. 1999); other suggestions are discussed by, 
e.g., Beck et al. (1996), Kulsrud (1999), Field et al. (1999) 
and Blackman & Field (1999). In particular, Blackman & 
Field (2000) argue that the i?m-dependent quenching seen 
in the simulations of Cattaneo & Hughes (1996) is a con- 
sequence of helicity conservation when using closed or pe- 
riodic boundaries, while simulations with open boundaries 
by Brandenburg & Donner (1997) (see also Brandenburg 
2000) do not show this effect. 
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The aim of this letter is to demonstrate that with open 
boundaries the scenario of Vainshtein and Cattaneo re- 
sults in basically the same estimate for the equilibrium 
magnetic field strength as is given by the naive viewpoint. 

The essence of our arguments can be presented as fol- 
lows. According to Vainshtein and Cattaneo, the suppres- 
sion of dynamo action by the small-scale magnetic field 
that is generated together with the large-scale is connected 
with the magnetic helicity of the small-scale magnetic 
field. Because the total magnetic helicity is an inviscid in- 
variant of motion, the magnetic helicity of the small-scale 
magnetic field can be connected with the magnetic helicity 
of the large-scale magnetic field. The governing equation 
for magnetic helicity has been proposed by Kleeorin and 
Ruzmaikin (1982; see the discussion by Zeldovich et al., 
1983), investigated by Kleeorin et al. (1995) for stellar 
dynamos, and self-consistently derived by Kleeorin and 
Rogachevskii (1999). During nonlinear stages of the dy- 
namo, the a-effect is thought to be determined by the 
hydrodynamic and magnetic helicities, so a closed system 
of equations can be obtained for the evolution of the mag- 
netic field and the a-coefficient (sec below, Sect. 2). This 
governing system (with helicity locally conserved) leads to 
magnetic field behaviour which is consistent with the pre- 
diction of Vainshtein and Cattaneo (we are grateful to M. 
Reshetnyak, who provided us with the relevant numerical 
results, which will be published elsewhere). 

We stress that Eq. ([I|) takes into account the local 
helicity balance at a given point inside the galactic disc 
\z\ < h, r < R, where r,<p,z are cylindrical coordinates. 
However, the kinematic galactic dynamo is impossible 
without a turbulent flux of magnetic field through the sur- 
face \z\ = h (see, e.g., Zeldovich et al. 1983, Ch. 11). It is 
more than natural to believe that this flux can transport 
magnetic helicity to the outside of the disc. The methods 
of Kleeorin and Rogachevskii (1999) allow us to introduce 
the corresponding term into the governing equations for 
the galactic dynamo. We demonstrate by numerical sim- 
ulations, and to some extent analytically that this term 
leads to a drastic change in the magnetic field evolution. 
Now the steady-state large-scale magnetic field strength 
is approximately in equipartition with the kinetic energy 
of the interstellar turbulence. 



2. Equations for magnetic helicity 

Following Kleeorin and Ruzmaikin (1982), we parameter- 
ize the back-reaction of dynamo generated magnetic field 
in terms of a differential equation for the a-coefficient, 
using arguments from the magnetic helicity conservation 
law. It is necessary to introduce the large-scale vector po- 
tential A, small-scale vector potential a, and the corre- 
sponding representations for the magnetic fields, B and 
b. We then write the total magnetic field as H = B + b, 
and the total vector potential as A — A + a, thus de- 
composing the fields into mean and fluctuating parts. The 
equation for the vector potential A follows from the in- 



duction equation for the total magnetic field H 

dA/dt = v x H - 7] curl H + Vip , (1) 

where v = V + u, and V = (v) is the mean fluid veloc- 
ity field, r\ is the magnetic diffusion due to the electrical 
conductivity of the fluid, ip is an arbitrary scalar function. 
Now we multiply the induction equation for the total mag- 
netic field H by a and Eq. ([!]) by b, add them and average 
over the ensemble of turbulent fields. This yields an equa- 
tion for the magnetic helicity x h — ( a ■ b) in the form 

dx h /dt + V • F = — 2(u x b) • B — 2?7(b • curl b) , (2) 

where F = (2/3) Vx h + (ax(u x B)) - ^(axcurlb) + 
(ax(u x b)) — (hip) is the flux of magnetic helicity. The 
electromotive force for isotropic and homogeneous turbu- 
lence is 

(u x b) = aB — ?7rcurlB , (3) 

where r\T is the turbulent magnetic diffusivity, and it is as- 
sumed that a is the total alpha-effect which at the nonlin- 
ear stage includes both the original hydrodynamical, and 
the magnetic, contributions. Note that the magnetic part 
of the a effect is proportional to the magnetic helicity, i.e. 
a h _ J Qftjrqrp (gee, e.g., Kleeorin and Rogachevskii, 
1999), where p is the density. The simplest form of the 
magnetic helicity flux for an isotropic turbulence is given 
by F = Yx h ', where V is the mean fluid velocity, e.g. that 
of the differential rotation (see Kleeorin and Ruzmaikin, 
1982; Kleeorin and Rogachevskii, 1999). Thus, the equa- 
tion for the magnetic part of the a effect in dimensionless 
form is given by 

^ + ^+V-(V X ' 1 ) = 4(V0 2 (^B ■ curlB-aS 2 ) , (4) 

(see Kleeorin and Ruzmaikin, 1982), where I m 100 pc is 
the scale of turbulent motions. We adopt here the stan- 
dard dimensionless form of the galactic dynamo equation 
from Ruzmaikin et al. 1988; in particular, the length is 
measured in units of the disc thickness h, the time is mea- 
sured in units of h 2 /rjx and B is measured in units of the 
equipartition energy B cq = ^/4npu. Here u is the char- 
acteristic turbulent velocity in the scale I, T]t — lu/3, 
T = (l/3)(l/h) 2 Rm and R a — la*/r]T, where a h and a 
are measured in units of a* (the maximum value of the 
hydrodynamic part of the a effect). For an axisymmetric 
dynamo V • (V% ) = 0. 

When da h /dt = and R~ X B ■ curlB < aB 2 , Eq. (@) 
yields a — a v /[l + (4/3)RmB 2 } (see, e.g., Vainshtein and 
Cattaneo , 1992). However, the latter equation is not valid 
for galaxies because da h /dt 3> a h /T. In addition, the 
condition -fi^B • curlB <C aB 2 seems not to be valid for 
galaxies. 

Equation (^) has been later reproduced, e.g. by 
Gruzinov and Diamond (1995). However, although this 
equation has never been included into detailed galactic dy- 
namo calculations, nevertheless its qualitative properties 
are more or less clear. Provided that dissipative losses are 
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taken into account, Eq. (||) leads to the same type of be- 
haviour as that obtained by the ad hoc prescription of the 
result of Vainshtein and Cattaneo (1992), i.e. the steady 
state strength of magnetic field is about B ca Rm~ x l 2 (see, 
e.g. Field et al., 1999). The real advantage of Eq. (^) is 
the fact that it is derived from first principles rather than 
prescribed ad hoc. If the dissipative losses in Eq. ([|) are 
neglected, the magnetic field decays for t — > oo. We stress 
that Eq. ^ contains a large factor i(h/l) 2 ~ 100 typi- 
cally. 

Kleeorin and Rogachevskii (1999) extended the calcu- 
lations to include a flux of magnetic helicity. Based on 
Eq. (13) of that paper, the approximate relation 



where 



<P(B) 



UB 2 



arctan(v / 8-B) 



2B 2 [l 



^ = 4^) [(B-curlBi?- 1 -^)^) 

+ ^(a v (z)4>{B)B 2 hh(z))] 



16B 2 + 1285 4 ln(l + (8B 2 )" 1 )] 



Thus 4>{B) = \/(AB 2 ) for B > 1/V§ and = 1 - 
(48/5)B 2 for B < 1/V8. The function is derived 
by Rogachevskii and Kleeorin (2000). Note that in a 
more simplified model of turbulence the function 4>{B) = 
M B ) = (3/8B 2 )(l - arctan(V8B)/V8 J B) (see Field et 
al. 1999). We stress that the qualitative behaviour of the 
model does not depend on these uncertainties in estimates 
for the scaling functions <f> and <j>\. 

Now we insert the a-coefficient given by Eq. (^J) into 
local disc dynamo problem to obtain the following equa- 
tions: 



(5) 



b" 



dBp 
dt 



Db r + B'l 



can be formulated. In Eq. (||), fi(z) describes the in- 
homogeneity of the turbulent diffusivity, and we define 
f(z) = a v (z) fi(z) . The profile f(z) depends on details of 
the galactic structure. Also, a(B) is the total a effect and 
a = a v (j){B) + a h (pi(B), where B = |B|. Here a v is the 
hydrodynamic part of the a effect, with a v 4>(B) its mod- 
ification due to nonlinear effects. Correspondingly, a h is 

the magnetic part of the a effect, and a h (j>i(B) is the mod- B • curlB = R a (B^b' r — b r B'i) 7 



(7) 

(8) 



M 
(0) 



(here B r = R a br). We can then obtain the steady-state so- 
lution of Eqf 
coordinates 



lution of Eqs. (uh and (H). Recognizing that in cylindrical 



ification caused by nonlinear effects (see Rogachevskii and 
Kleeorin, 2000). <j>i(B) = (3/8B 2 )(l-arctan(V8B)/V8S) 
and the function <fi(B) is defined below. The magnetic part 
of the a effect is proportional to the magnetic helicity, i.e., 
a h _ ^ /l /(i8 7r ^ T p) (see, e.g., Kleeorin and Rogachevskii, 
1999). For galaxies the term a h /T is very small and can 
be dropped. The gauge conditions V ■ A = V • a = have 
been used; our results can be shown to be gauge invariant 
(see Berger and Ruzmaikin, 2000). 

The last term in Eq. (JsJ) is related to the turbulent flux 
of magnetic helicity. This turbulent flux is proportional 
to the hydrodynamic part of the a effect and the turbu- 
lent diffusivity (see Kleeorin and Rogachevskii, 1999). The 
turbulent flux of magnetic helicity serves as an additional 
nonlinear source in the equation for the magnetic part of 
the a effect and it causes a drastic change in the dynamics 
of the large-scale magnetic field. 

For simplicity we replace the flux divergence in the 
right hand side of Eq. (||) by a decay term, i.e. we replace 
Jf- by 1/h (in principle, there is no problem in treating 
this point more carefully). 



3. The equilibrium magnetic field configuration 

We now present some asymptotic expansions for galactic 
dynamo models with the nonlinearity (||). First of all, we 
recognize that, because of the large parameter A{h/l) 2 in 
the right hand side of Eq. ([|) , we can take 



(9) 



we obtain for fields of quadrupole symmetry (cf. Kvasz et 
al., 1992) 



B'l' + Da(B)B 4 







(10) 



in a steady state. The corresponding equation in kinematic 
theory reads 

B'l' + Da Q B^ = Q. (11) 

Substituting (^) into ([!(]) we obtain, 

B'I'B 2 + DB 4> [f{z)4>{B)B 2 + R^B ■ curlB] = 0. (12) 

Using Eq. (^J) we rewrite Eq. ( |l2[ ) in the form 

B';\B 2 - B 2 ) + B^B't + Df(z)^(B)B 2 } = 0. (13) 

For the att dynamo B w B^. This assumption is justified 
if \D\ » R a , i.e. \R U \ » 1. Eq. (|l|) then becomes 

B"B' + Df(z)4>(B)B 2 = 0, (14) 



Note that Eq. ( |14| ) differs from Eq. flll]), arising from 
kinematic theory. For the specific choice of helicity pro- 
file f(z) — sin7rz and negative dynamo number D, there 
is an explicit steady solution, if we assume B 2 « B 2 ^ (re- 
member that also B 3> 1/V8, i-e. super-equipartition) , of 
the form 



2y/\D\ irz 



a(B) = f(z)cj)(B) + R^B^B ■ curlB, 



(G) 



Br — — 



y/nRa 



1TZ 

-.B eq cos — , 



(15) 
(16) 
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where we have restored the dimensional factor i? cq - (Note 
that B • curlB = for this approximate solution.) This 
solution is remarkably close to the results from the naive 
Ansatza = a (l—(B / B cq ) 2 ) ora = a / (l + (B/ B eq ) 2 ), or 
the model of Moss et al. (1999). For example, the pitch an- 
gle of the magnetic field lines is p = — arctan (7r 2 /4|i? w |) 
14° for \D\ = 10 and R a = 1. 



4. Numerical results 

We verified numerically that the initially weak magnetic 
field approaches the equilibrium configuration ( |l5| ) with 
accuracy 1% for \D\ > 1000, and an accuracy of 50% for 
|D| > 10. As is anticipated in the previous section, the 
equilibrium magnetic field near to the generation thresh- 
old value is more complicated. The threshold value for 
the nonlinear solution of Eqs. (Q) and (||) is D w —3.14, 
while the linear threshold value is D « —8. This is be- 
cause the nonlinear solution arranges itself so that the 
term B • curlB/5 2 , in a (see Eq. (||)) is of order 1. Thus, 
for the nonlinear solution with D — —8, the maximal value 
of a is about 1.25, whereas for D = —5, the maximal value 



is about 1.76. For |D| ~ 10 we obtain numerically 



B (O) ~ 0.23|L> - D c 



I).") 2 



(17) 



where D cl - is the nonlinear threshold value. As \D — D cr \ 
increases towards 10, the slope increases slightly, but Eq. 
( |l7j ) remains a reasonable estimate. (Note that accurately 
estimating the exponent in Eq. (y_7|), and subsequently, 
is a quite delicate matter even in this one-dimensional 
problem, and that the quoted figures may be uncertain in 
the last digit.) 

This result is robust under variations of the helicity 
profile. For f(z) = z we get in the nonlinear case D CI = 
—7.49, while the linear threshold value is D CI = —12.5 
and B (O) « 0.15|D - L> C r| a5 ° near D = D cr , i.e. again a 
square root dependence to within the errors of our proce- 
dure. Further, with f(z) — z/\z\ in \z\ > 0.2, and a smooth 
interpolation to zero in \z\ < 0.2, we find D CI w —2.41 
and £^(0) w 0.25|D — -D cr | a50 , again closely the same 
dependence. (In this case the linear threshold value is 
D cr = -6.53.) 



5. Discussion 

We have demonstrated that the nonlinear evolution of the 
helicity following from Eq. (^) gives a basically different 
type of galactic magnetic field evolution to that follow- 
ing from Eq. (||). Eq. (^), being based on local helicity 
conservation, results in magnetic field decay, after a stage 
of kinematic growth. If the molecular diffusivity of the 
magnetic field is taken into account, this decay is followed 
by a stabilization at a very low magnetic field strength, 
corresponding to the estimate of Vainshtcin and Cattaneo 
(1992). The scenario of magnetic field and helicity dynam- 
ics can then be described as follows. Large-scale dynamo 
action produces large-scale magnetic helicity. Due to the 



local conservation of helicity, suppression of field genera- 
tion results. An equilibrium is possible if molecular diffu- 
sivity is present, so the equilibrium magnetic field strength 
is very low. 

Equation (|J) allows for the transport of helicity, so 
the local value of the helicity changes during magnetic 
field evolution. The scenario of magnetic field and helic- 
ity dynamics can be presented as follows. As usual, mag- 
netic helicity of the large-scale magnetic field is produced, 
however the total magnetic helicity is not now conserved 
locally, but the magnetic helicity of the small-scale mag- 
netic field is redistributed by a helicity flux. The equilib- 
rium state is given by a balance between helicity produc- 
tion and transport. The helicity conservation law now ex- 
presses the conservation of an integral of the helicity over 
the galactic disc. However this conservation law is triv- 
ial, because the integral vanishes identically as helicity is 
an odd function with respect to z. Now the equilibrium 
strength of the large-scale magnetic field is of order that 
of the equipartition field: this is our main result. 
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